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ABSTRACT 



In situ measurements of salinity and temperature fluctu- 
ations at depths to 14 meters indicate distinct dependences 
at different times of the day. The variance of the salinity 
fluctuations decreased with increasing depth, but was greater 
just after sunrise than just prior to sunset. The variance 
of the temperature ■ fluctuations decreased with increasing 
depth just prior to sunset, but increased with depth immedi- 
ately after sunrise. The correlation length of the sound 
index of refraction was calculated by using the variance of 
the sound velocity fluctuations, and the variance of sound 
amplitude modulation in the theory of Mintzer. This analysis 
shows that microstructure patch size increases approximately 
linearly with depth. The power spectral densities of the 
salinity, temperature and sound velocity fluctuations show 
peaks of energy corresponding to dominant ocean wave 
frequencies . 
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I. INTRODUCTION 



A. HISTORY 

The study of sound propagation in the ocean and the ef- 
fects of oceanographic parameters (particularly temperature) 
on it has been the subject of considerable research since 
early World War II. Wilson [Ref. 1] formulated an equation 
describing the relationships between sound velocity, tem- 
perature and pressure in distilled water, and then extended 
his work to an oceanic regime by utilizing filtered seawater 
diluted with distilled water to obtain desired salinity 
dependence [Ref. 2]. Since Wilson's efforts in this field, 
many investigators have attempted to revise the sound speed 
equation with newer techniques and theory. The reasons for 
revised formulations included, among others, the use of new 
type velocimeters [Ref. 3] the suspicion that standard sea 
water did not maintain a standard salinity to sound speed 
relationship [Ref. 4] , the desire to obtain simpler expres- 
sions, and the use of a more realistic range of values for 
the oceanographic parameters [Ref. 5] . 

While these investigations into the parametrical rela- 
tionships to sound velocity were proceeding, other investi- 
gators attempted to measure the magnitude of the fluctuations 
of the important parameters affecting sound propagation in 
seawater. Liebermann [Ref. 6] studied the effects of tem- 
perature inhomogeneities on sound propagation. He found that 
thermal inhomogeneities do have an effect on the scattering. 
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intensity, and refraction of the sound transmission. Shon- 
ting [Ref. 7] also measured the thermal microstructure, his 
results indicating the importance and magnitudes of tempera- 
ture fluctuations in the ocean environment. 

B. SIGNIFICANCE OF PROBLEM 

The characteristics of sound propagation in a deep ocean 
environment (defined in this paper as the depth below the 
thermocline) are well known and predictable. This is due 
to the fact that in the deep water regime the fluctuations 
of temperature and salinity are minor. As a result of this 
predictability, sound is used properly and constructively 
in deep waters. For the United States Navy this specifically 
relates to optimum usage of sound in the field of anti- 
submarine warfare . 

However the propagation of sound in the shallow ocean 
environment and in the near surface region is not as pre- 
dictable, and as a result seriously hinders the user and 
the planner in their efforts to utilize sound to its best 
advantage. The near surface region usually plays a part in 
long range propagation of sound. 

Sound velocity amplitude and phase in the region between 
sound source and sound receiver change in the near surface 
regions. They change due to small scale variations in oceano- 
graphic parameters including possibly the presence of micro- 
bubbles and organic content in the regime. 
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C. OBJECTIVES OF THIS RESEARCH 

The main motivation in this research was to statistically 
describe the oceanic parameters of temperature and salinity, 
and to relate their statistics to the statistics of sound 
velocity in a near surface environment. It was hoped to ob- 
tain enough data to enable one to draw conclusions concern- 
ing the variation of the fluctuations with respect to depth 
and time of day. If such concrete conclusions could be 
reached in the near surface region, then the users of sound 
either in oceanographic research or in present anti-submarine 
warfare would be afforded new tools with which to update the 
state-of-the-art in sound use and understanding. 



13 



IX . EXPERIMENTAL PROCEDURE 



A. OCEANOGRAPHIC ENVIRONMENT 

1. Location 

Due to its advantageous location of being rigidly 
fixed in a shallow water environment the N.U.C. Oceano- 
graphic Research Tower, approximately one mile off Mission 
Beach, California, was chosen as the location for the pro- 
ject. It allows for continuous data accumulation and 
monitoring as well as some analysis in a near laboratory 
controlled environment. 

Figure 1 is an illustration of the tower. It is 
fixed by supporting pins driven 63 feet into the ocean 
floor. Electrical power is supplied from shore thereby 
guaranteeing stable voltage and frequency. The tower is 
located in 60 feet (18 m) of water with a sandy bottom 
and is free of water traffic most of the time. The upper 
level of the two level tower was the location for elec- 
tronic support equipment. The lower was for handling 
equipment. 

2 . Time of Year 

The experiment was conducted from the afternoon of 
21 October 1971 through the morning of 22 October 1971. 

3. Meteorological, Oceanographic, and Biological 

Conditions at Time of Project 

a. Weather State: clear to light haze 

b. Wind: 21 October - 300-320 from 7 to 9 kts 

22 October - 065-110 from 2 to 8 kts 
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FIGURE 1. NUC Oceanographic Research Tower. 
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c. Wind Waves: 0 to h ft from crest to trough 

d. Swell: 2 ft from Northwest 

e. Bubbles: 1 square foot in vicinity of tower 

stanchions during all runs. During run #4 bubbles were 
present on the surface in patches of size 20 ft by 2 ft 
near vicinity of tower. 

f. Biologies 

(1) Seaweed . Small amount of sea grass on 
surface during run #4. 

(2) Animal' Life . Seal in vicinity during run 
#6, porpoise and seals in vicinity during runs 7 and 8. 

g. Temperature Structure: thermocline varied from 

depths between 40 ft to 60 ft below the surface (see 

Fig. 2) 

h. Salinity: Nansen Cast results, as determined 

from Hytech model 621 salinometer, found in Table III. 

B. EQUIPMENT USED DURING EXPERIMENT 
1. Bissett-Berman STD Model 9006 

a. Temperature Sensor 

(1) Time Constant: 0.35 sec 

(2) Output: 0-10 mv DC 

(3) Temperature Range Used: 14 - 19° C 

(4) Accuracy: ±0.02° C 

b. Salinity Sensor 

(1) Time Constant: 0.35 sec 

(2) Output: 0-10 mv DC 
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Depth in Feet 



Temperature (°F) Temperature (°F) 





0800 22 Oct 1971 



0900 22 Oct 1971 



FIGURE 2. Bathythermograph Traces. 
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(3) 


Salinity Range Used: 


37.5 - 39 .'5 ppt 
(before correction) 


(4) 


Accuracy: ±0.03 ppt 





c. Detailed Description Found in Appendix A 

2. Thermistors ' (3 ) 

a. Time Constant: 150 milliseconds 

b. Output: .0553 volts per °C 

c. Temperature Range: Thermistors nulled using 

Wheatstone Bridge Circuit prior to start of each run. 

d. For further description see NPS thesis of LCDR 
Duchock, March 1972. 

3. Ramsay Corporation Mark I SVTD 

a. Time Constant: 160 microseconds 

b. Output: 0 - 10 v DC 

c. Sound Velocity Range: 1400 - 1600 m/sec 

d. Accuracy: ±.01 m/sec 

e. For further description see NPS thesis of LCDR 
Duchock, March 1972. 

4 . Wave Height Sensors 

a. Baylor Wave Profile Recorder System 

(1) Output: 50 millivolts per foot of height 

(2) Accuracy: 1% or reading 

(3) For further description see NPS thesis of 
LT Bordy, March 1972. 

b. Interstate Electronics Corporation Pressure Wave 
Gauge 

(1) Output: .25 v/psi 
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(2) Accuracy: 2% of reading 



5. Particle Velocity Sensor 

a. Engineering Physics Company Electromagnetic 
Current Meter 

(1) Output: variable voltage per meter/second 

(2) Accuracy: 1% of full scale reading 

b. For further description see NPS thesis of LT 
Bordy, March 1972. 

6 . High Pass Filter #1 - Krohn-Hite Model 3340 

a. Gain: OdB 

b. High Pass Filter Cutoff: 0.01 HZ 

7. Band Pass Filter #2 - Krohn-Hite Model 330-A 

a. Gain: OdB 

b. Low Pass Cutoff: 0.02 HZ 

High Pass Cutoff: 2000 HZ 

8. Band Pass Filter #3 - Velocimeter 

a. Gain: OdB 

b. Low Pass Cutoff: 0.02 HZ 

High Pass Cutoff: 2000 HZ 

9 . Amplifier #1: ' Preston Scientific Model 8300 xwb 

a. Gain: 500 

10 . Amplifier #2: Preston Scientific Model 8300 xwb 

a. Gain: 1 

11 . Amplifier #3: ' Hewlett-Packard Model 2470 

a. Gain: 30 

12 . Amplifier #4: Hewlett-Packard Model 2470 

a. Gain: 100 
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13. Amplifier #5; Hewlett-Packard Model 2470 

a. Gain: 100 

14 . Amplifier #6': ' Hewlett-Packard Model 24 70 

a. Gain: varied 100-300 

15 . Sangamo Magnetic Tape Recorder Model 3500 

a. • 14 channel with 11 channels allocated for 
oceanographic parameter data were utilized. 

b. Channel Allocation: 



(1) 


Ramsay Temperature Sensor (not 


used) 


(2) 


Ramsay Velocimeter 




(3) 


Bissett-Berman Salinometer 




(4) 


Turbulent Velocity (horizontal 


component) 


(5) 


Turbulent Velocity (vertical component) 


(6) 


Baylor Wave Height 




(T) 


Sound Amplitude Modulation 




(8) 


Pressure Wave Height 




(9) 


Thermistor #1 




(10) 


Bissett-Berman Temperature (not 


. used) 


(11) 


Thermistor #2 




(12) 


Thermistor #3 




(13) 


Sound Phase Modulation 




(14) 


Voice (not used) 





C. EXPERIMENTAL CONFIGURATION 
1. Schematics 

a. Bissett-Berman STD - See Fig. 3 

b. Ramsay Velocimeter - See Fig. 4 

c. Temperature Thermistors - See Fig. 5 
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FIGURE 3. Bissett-Berman STD Flow Diagram. 
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FIGURE 4. Ramsay Sound Velocimeter Flow Diagram. 
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FIGURE 5. Thermistor Flow Diagram. 
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2. Geometry of Sensors 

The sensors' geometric relationships were deemed to 
be of major importance in view of potential acoustical and 
physical interferences which could have had an adverse ef- 
fect on the experimental results. In addition to the attempt 
to describe the oceanic regime in terms of the oceanic 
parameters two acousticians simultaneously conducted ex- 
periments related to sound amplitude modulation, sound 
dispersion, and phase fluctuations. (See NPS thesis of 
LCDR W. J. Smith, Jr., December 1971 and thesis of LCDR 
Juergen Rautmann, December 1971.) 

In order to reduce the effects of mutual interfer- 
ence between the sensors especially in the path of sound 
propagation between hydrophone and receiver, a physical 
separation of the sensors was mandatory. The geometrical 
separations between the sensors as shown in Figs. 6, 7 were 
used. Photographs of the Bissett-Berman STD and the geo- 
metrical relationships of the other sensors to it appear 
in Figs. 8-11. 
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FIGURE 6. Geometrical Relationships Between Sensors. 
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FIGURE 7. Geometrical Separations Between 

Oceanographic Sensors and Acoustic 
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III. DATA COLLECTION AND RUN DESCRIPTION 



Data collection and run description will be found on page 56 under 



Appendix B. 
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IV. DATA REDUCTION 



A. INTRODUCTION 

The analog data were recorded for each run on the 
Sangamo Model 3500 magnetic tape recorder operated at a 
speed of 1 7/8 ips using FM electronics . The recorded 
data were then transferred to a Brush Mark 200 Strip-Chart 
recorder. In some instances, in order to produce better 
resolution, the recorded data were transferred to a Brush 
Mark 2 Strip-Chart recorder through a Krohn— Hite Model 
3340 filter set on low pass-max flat (see Fig. 12) . The 
printed analog data were converted to digital data on 
magnetic tape using the Fleet Numerical Weather Central, 
Point Pinos, California tracing digitizer. The taped data 
were transferred to punched IBM data cards compatible for 
statistical analysis on an IBM 360-67 digital computer. 

B. MAGNETIC TAPE TO STRIP-CHART RECORDING 

The recorder playback speed was seven and one half 
inches per second, which is an increase of four times over 
the data collection speed. The strip-chart recorder was 
run at a speed of five millimeters per second. The speed 
differentials between the record, playback, and strip-chart 
phases were programmed in the IBM program used for data 
analysis. The sensitivity setting on the strip-chart re- 
corder was important for data presentation and analysis . 

The sensitivity settings varied from run to run and with 
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FIGURE 12. Representative Oceanographic Fluctuations. 
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each parameter. The settings used can be found in 
Appendix B . 

C. ANALOG TO DIGITAL TAPE 

In order to utilize the IBM 360 computer to determine 
the statistics of the fluctuations , the analog data on the 
strip-chart record had to be converted to digital form. In 
order to effect the conversion use was made of the Calma 
tracing digitizer at Fleet Numerical Weather Central's 
Point Pinos, California facility. The tracing digitizer 
transfers analog data to digital data on magnetic tape at 
a sampling rate of 100 samples per inch in both the x and 
y directions . The length of each run analyzed amounted to 
approximately the first nineteen minutes of the run while 
the number of data points varied depending upon the posi- 
tioning of the strip chart on the digitizer. It was im- 
portant for future correlations that the tracing of all 
parameters began at the same starting point. This was 
ensured by starting at the sharp line transient signaling 
the start of a particular run. 

D. CONVERSION OF TAPED DIGITAL DATA TO DATA CARDS 

Conversion of taped digital data to a form compatible 

for use on the IBM 360 Digital Computer, i.e. data card 
form, was accomplished using Fleet Numerical Weather 
Central's Computer Facility. Fourteen data points were 
recorded on each IBM card. The number of data points per 
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run varied depending on the positioning of the strip chart 
on the Calma Plotter at Point Pinos, but the average run 
produced approximately 5400 data points. 
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V. DATA ANALYSIS 



A. INTRODUCTION 

Attention was focused entirely on analysis of the 
statistics of the fluctuating parameters. Two arguments 
led to this method of analysis: first of all, the statis- 

tics of the oceanographic variables are important in deter- 
mining the statistics of sound propagation. Since the 
major objective of the experiment was to define the medium 
in order to determine its effect upon the propagation of 
sound the statistical method of analysis was the natural 
approach. The second reason for using statistics as the 
major tool for analysis was due to the relatively wide 
separation between sensors which was required to prevent 
instrumental interactions. In order to enable one to dis- 
cuss the interrelationships between spaced sensors it was 
necessary either to correct for the lack of idential point 
location or to consider the statistical relations, assum- 
ing that the medium was statistically homogeneous in the 
region of the sensors. The latter is the easier choice. 

B. STATISTICAL PROGRAMMING 

The programs used in the data analysis appear in 
Appendix C . 

Program A analyzed the data for mean, variance, normal- 
ized autocorrelation function, and power spectral density. 
The standard Blackman- Tukey method was used to compute 
power spectral density with a Parzen window as the spectral 
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window. Program A was modified to compute variance by two 
methods. The first method calculated variance by integra- 
ing the power spectral density. A subroutine (aver) was 
then added to calculate variance directly from the time 
series data. 

Program B was used to output data from which a histogram 
could be plotted of the number of times a particular value 
of a parameter occurred versus the value of that parameter. 

In utilizing program A several important inputs were 
necessary : 

NTS: number of data points 

MLAG: number of time lags = 10% (NTS) 

DT : time increment in seconds/sample. This parameter 

took into account the different driving speeds of the strip- 
chart recorder, playback and record modes of the tape re- 
corder, and the analog-to-digital digitizer. 



DT =- 



inch 100 samples 

(5/4 mm/sec) ( ) ( ) 



25 . 4mm 
= 0.2032 sec/sample 
Af: frequency increment = 



inch 



2 MLAG (DT) 



FBHZ : lowest frequency 

FEHZ : highest frequency = Nyquist Frequency 

1 



FEHZ = 



= 2.46 HZ 



2 ( .2032) 

Calibration factor: variable on each run 
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The raw data on the IBM cards had the units of inches. 
The calibration factor converted inches to volts per milli- 
meter by incorporating the Brush recorder sensitivity set- 
ting into the program input thereby giving the data units 
of volts per inch. When each data point was multiplied by 
the calibration factor the results produced were in voltage 
units . 

# volts 25.4 millimeters volts 

C.F. = ( ) = "X" 

millimeter inch inch 

C. STATISTICAL ANALYSIS 

1 . intro due ti on 

In order to convert the program output in volts to 
the correct oceanographic units and values, conversion 
factors for sensor voltage sensitivities, amplification 
factors and filter gains were used. On location monitoring 
of the data combined with post-experimental review in- 
dicated some runs higher in quality and accuracy than 
others. For such reasons runs 3,4,5,6,7,9,10 were analyzed 
for salinity and temperature, and runs 4, 5, 6, 8 for sound 
velocity. The conversion factors used were: 
thermistors: 5.53 v/°C 

salinity: 2.5v/%. 

m 

sound velocity: 1.5 v/ 

sec 

2 . Variances of Oceanographic Fluctuations 

The values of the variances of the oceanographic 
parameters as computed by integrating the power spectral 
densities are indicated in Table I. 
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The values of the variances of the oceanographic 
parameters as computed directly from the time series data 
are listed in Table II. 

It was noted that there was only a very slight dif- 
ference between the values of variance computed by both 
methods. However the variances used in subsequent analysis 
were those computed directly from the data thereby ensuring 
a higher degree of accuracy. Additionally all subsequent 
analysis' was performed using temperature data measured by 
thermistor #1 since the quality of the fluctuating signal 
from that sensor was the best and most consistent throughout 
all runs. (Thermistor #2 data used during run 4 since 
thermistor 2 was operating reliably at that time) . Figures 
13-15 are graphical displays of variances versus time of day 
with depths at which runs were made noted in parentheses. 

A survey of the literature indicates that this ex- 
periment was unique in its attempt to establish statistical 
relationships between the described oceanographic parameters. 
As a result little direct comparison with previous results 
is possible. 

3. Variance of Temperature Fluctuations 

Two past investigators have studied temperature fluc- 
tuations in the near surface region on approximately the same 
time scale as was used in this experiment. Shonting used a 
double-beaded thermistor manufactured by Fenwall Electronics, 
Inc., Framingham, Massachusetts, to observe thermal micro- 
structure. He conducted his measurements in the Bahamas in 
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Run 


Time 


Depth 


Sound Vel. 
Variance 


Temp # 1 
Variance 


Temp # 2 
Variance 


Salinity 

Variance 


3 


21 OCT 
1616 


9 « 0m 




.0187 (°C) 2 


. 0121 (°C) 2 


. 0016 (%.) 2 


4 


1648 


13.9m 


. 0294 (M/S) 2 


. 0167 


_ 


. 0013 


5 


1728 


6 . 9m 


. 0003 


. 0230 


. 0002 


.0023 


6 


22 OCT 
0354 


4. 3m 


. 0072 


.0052 




. 0044 


8 


0650 


9 . 5m 


.0023 


.0013 




.0110 


9 


0725 


14 . 6m 




. 0076 




.0065 


10 


0802 


7 . 8m 




. 0017 


. 


. 0002 



TABLE I 



Run 


Time 


Depth 


Sound Vel. 
Variance 


Temp #1 
Variance 


Temp # 2 
Variance 


Salinity 

Varianae_ 


3 


21 OCT 
1616 


9 . 0m 




,0143 (°C) 2 


. 0093 (°C) 2 


.0016 (%•) 2 


4 


1648 


13.9m 


. 0293 (M/S) 2 ! 


. 0129 


. 


.0012 


5 


1728 


6 . 9m 


.0003 ! 


. 0184 


. 0002 


. 0023 


6 


\ 

22 OCT 

0354 


4.3m 


. 0072 


. 0042 


i ” 


. 0041 


8 


0650 


9 . 5m 


. 0021 


. 0010 




.0111 


9 


0725 


14.6m 


j 


. 0059 




. 0059 


10 


0802 


7 . 8m 




. 0013 


m , 


.0002 



TABLE II 
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August 1967. Analysis of his measurements taken between 
1323-1400 indicated representative values of variance of 
approximately 25 x 10“ 6 (°C) 2 at two meters depth and ap- 
proximately 72 x 10“ 8 (°C) 2 at 20 meters, a decrease of two 
orders of magnitude with depth down to 20 meters . Although 
the variance was much greater in the present experiment 
than in Shonting's, there was in the late afternoon period 
a one order of magnitude decrease in temperature variance 
with depth down to approximately 15 meters. 

Liebermann also recorded temperature fluctuations in 
the near surface region at a depth of 50 meters. A rough 
calculation of variance from Liebermann' s Fig. 1, which was 
given as a representative temperature record for July and 
August (time of day unreported) for the continental waters 
from Southern California to Alaska, resulted in a variance 
of approximately 17 x 10 -4 (°C) 2 . This figure is somewhat 
closer to data of Fig. 13 but the lack of time of day in- 
formation and depth dependence prevents conclusions from 
being drawn. 

Due to the maximum temperature gradient being effected 
just prior to sunset the strongest fluctuations of tempera- 
ture were noticed at all depths at that time. The variance 
decreased with increasing depth. However in the short 
period just after sunrise the temperature fluctuations were 
noticed to increase with depth. Sagar [Ref. 8] quotes, from 
a personal communication, Deacon, of the National Institute 
of Ocoar. . v. • same s: 
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Variance in (°C) 



FIGURE 13. 



Plot of Temperature 
Variance Versus Time 
of Day (Depth in 
Parentheses) . 
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Variance in (°/oo) 



FIGURE 14 



Plot of Salinity 
Variance Versus Time 
of Day (Depth in 
Parentheses) . 
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21 Oct 
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22 Oct 

Sunrise 0703 



Time of Day by Hour 
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Variance in (m/sec) 



FIGURE 15. 



CM 




Plot of Sound Velocity 
Variance Versus Time 
of Day (Depth in 
Parentheses ) . 
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This behavior could possibly have been caused by the "con- 
vective lag" phenomenon. James [Ref. 9] describes this 
behavior as occurring in the early morning just after sun- 
rise. Insolation increases to the point where it is equal 
to the outgoing heat loss from the water. However heat 
loss occurs exclusively at the surface, with consequential 
continued cooling of the surface. As insolation increases 
in the hours just after sunrise, surface cooling continues 
while the lower levels are slowly heated. Eventually, of 
course, the lag is overcome and surface heating begins 
causing growth of small cells with large fluctuations near 
the surface and large cells of small fluctuations at greater 
depths, as described by Deacon. 

4 . Variance of Salinity Fluctuations 

The analysis of the plots of salinity variance versus 
time showed a decrease in salinity fluctuations with depth 
during the late afternoon period. The fluctuations were 
smaller during the period just prior to sunset than just 
after sunrise (assuming that the 7.8m reading is in error) . 
The decrease in variance of salinity with increasing depth 
in the late afternoon was expected since a stronger salinity 
gradient would exist in the surface layer where evaporation 
was initiated. However, it is speculated that convective 
mixing caused a more nearly isohaline condition to exist in 
the water column, by sunset than existed just after sunrise 
thereby accounting for the difference in magnitudes of 
variances at those two times of the day. 
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In the section "Parametric Interrelations Between 



Salinity, Temperature, and Sound Velocity" the theoretical 
calculation of sound velocity variance is discussed, and 
further analysis of the difference between theoretical and 
experimental results is made. 

5. Probability Density Function Results 

Salinity and sound velocity information were put into 
Program B in order to establish the type of distribution of 
the populations sampled. In all cases the distributions were 
Gaussian in nature. Figures 16 and 17 are examples of the 
distributions . 



6 . Normalized Autocorrelation Function 

The normalized autocorrelation functions of the oceano- 
graphic parameters (see Appendix D) indicated a long correla- 
tion time for temperature data thereby indicating relatively 
well mixed thermal structure. Conversely the relatively low 
degree of correlation for both salinity and sound velocity 
indicated an inhomogeneous medium with respect to those two 
parameters . The peaks and troughs found on the salinity and 
sound velocity autocorrelation plots indicated the presence 
of patches of higher and lower salinity water as well as 
patchy areas of water yielding higher and lower values of 
sound velocity. 

7. Power 'Spectral Density 

Plots were made of power spectral densities of salin- 
ity, temperature, and sound velocity fluctuations (see 
Appendix E) . 
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FIGURE 16 
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Salinity Variance in (volts) 



o 

VO 





<N 



# 



46 



Sound Velocity Variance in (Volts) 



The nature of the plots indicated peaks of "energy" 
centered about frequencies of 0.17 and 0.08 Hertz. This 
coincided with the dominant wave frequencies existing at 
the time of the experiment. 
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VI . PARAMETRIC INTERRELATIONS BETWEEN SALINITY, TEMPERATURE, 



AND SOUND VELOCITY 

A. INTRODUCTION 

With hopes of shedding light on possible heretofore un- 
known variables capable of affecting sound velocity a com- 
parison was made between sound velocity variance derived 
from Wilson’s sound speed equation and that measured experi- 
mentally with the Ramsay Velocimeter. 

B. DERIVATION OF THEORETICAL SOUND VELOCITY VARIANCE FROM 
WILSON'S EQUATION 

1. Wilson's Equation 

C (t) = 1449 .2 + 4.623 T (t) - .0546 T 2 (t) 

+ 1.391 (S (t) - 35) 

2. Computation of Variance 

var C = 0 + (4.623) 2 var T - (.0546) 2 var (T 2 ) 

+ (1.391) 2 var (S) 

3. Dropping the Term, var (T 2 ), as Negligible Yields 
the Final' Result 

var C = 21.372 var T + 1.935 var S 

The utilization of the above formula yielded the 
following results, with the velocimeter results shown for 
comparison : 
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i 



RUN 



VARIANCE C 



VARIANCE C. 



3 

4 

5 

6 
8 
9 

10 



theory 
.3095 (m/sec.) 2 

.2770 (m/sec.) 2 
.3979 (xn/sec) 2 
.0968 (m/sec) 2 
.0437 (m/sec) 2 
.1371 (m/sec) 2 
.0273 (m/sec) 2 



velocimeter 

.0293 (m/sec.) 2 
.0003 (m/sec.) 2 
.0072 (m/sec) 2 
.0021 (m/sec) 2 



It was hoped that the results would indicate a theo- 
retical variance slightly smaller than experimental variance 
thereby implying the presence of other parameters such as 
bubble presence to which the 3 MHz acoustic velocimeter 
would have been insensitive. However, as seen above, the 
experimental variance was much smaller than that calculated 
from the variances of temperature and salinity. Analysis 
of the wide discrepency led to the conclusion that since 
all sensors appeared to have operated in an optimum manner, 
and salinity and temperature signals were strong and yielded 
good results, the sound velocity signal was suspect. The 
signal had been filtered and amplified and conceivably could 
have been altered, or else a mistake made in the data re- 
duction settings . In order to determine which values of the 
variance of sound velocity should be used for the remainder 
of the analysis a final check of the results of both methods 
was conducted. Part of the experiment involved the measure- 
ment of sound amplitude modulation (see NPS thesis of LCDR 
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W. J. Smith, Jr., December 1971). Stone and Mintzer [Ref. 
10] define the coefficient of variation (CV) as: 



CV 




|P| 



2 



where AP is the variation of the acoustic pressure relative 

to the mean value P. In his experiment LCDR Smith measured 

cr which he shows is related to CV by: 
s J 



(CV) 2 = 40 a s 2 

Chernov [Ref. 11] formulated the relationship: 



(CV) 2 = VT y 2 K 0 2 aL 
where , 

2tt 

K 0 = — 

X 

X = acoustical wave length 
y 2 = variance of sound velocity 
a = correlation length (Gaussian correlation function 
assumed) 

L = distance between transducer and receiver 
This relationship combined with the measured values of CV 
and sound velocity variance, allowed for the calculation of 
the correlation length: 



RUN 


DEPTH 

(m) 


a <from Tm) Ctheory) 


a (frOI ”(c“ ° expt> 


5 


6.9 


175.0 


2.699 x 10 5 


3 


9.3 


157.5 


- 


4 


13.9 


452.5 


4279 .9 
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The values of "a" calculated from velocimeter data are un- 
reasonably large. For this reason the theoretical sound 
velocity variances are used throughout the remainder of the 
analysis. A graphical plot of theoretical sound velocity 
variances versus time of day with depths in parentheses are 
shown in Fig. 15. 

It is interesting to observe that Skudrzyk [Ref. 12] 
also found patch sizes, which increased approximately 
linearly with depth in approximatley the same manner as 
shown above. However, Skudrzyk' s patch sizes were much 
larger than reported here, probably due to his acoustic 
range being much greater than in this experiment. 
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VII. CONCLUSIONS 



A. Although the near surface regime is one of complex and 
high variability, several conclusions can be made: 

1. Temperature fluctuations decrease with depth during 
the period just prior to sunset. 

2. Temperature fluctuations increase with depth just 
after sunrise, but are smaller in magnitude than the 
fluctuations just prior to sunset. 

3. Salinity fluctuations decrease with depth but are 
larger in magnitude in the period after sunrise. 

4. Sound velocity fluctuations normally follow the 
trend of temperature variations. However, the im- 
portance of salinity fluctuations cannot be neglected. 
In one instance the salinity fluctuations had as much 
an effect on sound velocity fluctuations as did tem- 
perature fluctuations . 

5. The maximum energy spectra of salinity and temper- 
ature fluctuations occurred at frequencies that agreed 
with dominant surface wave frequencies. 

6. The sizes of patches of sound velocity fluctuations 
were calculated from the statistics of sound amplitude, 
temperature and salinity fluctuations . The implied 
patch sizes increased approximately linearly with in- 
creasing depth. 
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VIII. 'RECOMMENDATIONS FOR FURTHER STUDY 



A. While the medium, was described in terms of temperature 
and salinity, a major question remains, "Are there other 
variables that could affect the propagation of sound in 
sea water?" In order to best achieve an answer to this 
question one should ensure that a correct experimental 
value of sound velocity is measured for comparison with 
sound velocity computed theoretically. LCDR Smith's and 
Rautmann's thesis work indicated the presence of bubbles 
in the medium that could have a marked affect on sound 
velocity . 

B. In order to complete the description of the medium, an 
analysis of the size of the "patches" of temperature and 
salinity should be performed by use of arrays of sensors . 
Application of knowledge of the drift velocity, in addition 
to the procedures used in the present experiment, would 
also enable patch size determination. 
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APPENDIX A 



GENERAL DESCRIPTION OF BISSETT-BERMAN MODEL 9006 SALINITY, 
TEMPERATURE AND DEPTH MEASURING SYSTEM [Ref. 13] 

The model 9006 STD is comprised of underwater sensors 
and deck equipment whose functions are to furnish power to 
the sensors and process data received from same. 

In situ measurement of salinity is determined by measur 
ing conductivity, temperature and pressure. An inductively 
coupled sensor enables detection of conductivity, which is 
also compensated for temperature and pressure effects, 
thereby producing an output totally dependent upon salinity 
The output shifts the frequency of a PARALOC signal pro- 
viding an FM analog of salinity. 

The temperature sensor uses a platinum resistance 
thermometer and a PARALOC which provides an FM analog of 
temperature . 

The depth sensor utilized a pressure transducer into 
which is incorporated a strain-gauge bridge circuit. The 
changes in resistance are converted to a frequency analog 
in the PARALOC. 

The signal mixer receives and regulates power from the 
deck equipment and sends it onto the underwater sensors . 

It then multiplexes and amplified the FM signals from the 
sensor and sends them to the deck equipment. 

The distribution amplifier amplifies the FM signal 
from the underwater units, and separates data signals from 
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the DC power which enters the distribution amplifier from 
the power supply rack. 

Each underwater sensor has a separate discriminator 
which separates the mixed data signal into an individual 
channel for recording availability. The discriminator 
filters the desired signal and rejects others. The accepted 
signal is amplified, squared, and converted to a 10 millivolt 
DC level which is the output of the recorder. Dial readouts 
enable continual monitoring of this signal. 

The Power Supply is a model 8600 and provides a constant 
voltage of 26.5 volt DC to the deck units, and a constant 
current of 150 milliamps at a maximum of 140 volts DC to the 
underwater units. The underwater rack is stainless steel. 

The sensors and mixer are mounted in individual shock- 
absorbent brackets . 

A single conductor cable provides for physical connec- 
tion between sensors and deck equipment. 
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APPENDIX B 



DATA CUMULATION AND RUN DESCRIPTION 

The output DC voltage from the oceanographic sensors 
was recorded on the Sangamo magnetic tape recorder as 
explained in "Experimental Configuration" section. The 
runs during which data were collected are described below 



RUN 


DEPTH 


DATE 


START TIME 


STOP TI] 


1 


4 . 3m 


21 Oct 


1419 


1429 


2 


4.2m 


21 Oct 


1530 


1550 


3 


9 ,0m 


21 Oct 


1616 


1636 


4 


13.9m 


21 Oct 


1648 


1708 


5 


6.9m 


21 Oct 


1728 


1748 


6 


4 . 3m 


22 Oct 


0354 


0448 


7 


4 . 3m 


22 Oct 


0546 


0622 


8 


9 .5m 


22 Oct 


0650 


0715 


9 


14 . 6m 


22 Oct 


0725 


0749 


10 


7 . 8m 


22 Oct 


0802 


0821 


11 


8 . 2m 


22 Oct 


0832 


0852 


12 


5 . 7m 


22 Oct 


0929 


0945 



Run 3 



salinity 
thermistor 1 
thermistor 2 



Run 4 



SENSITIVITY SETTINGS 



50 millivolts/millimeter 
100 mi Hi volts /mil lime ter 
100 millivolts/millimeter 



salinity 50 millivolts/millimeter 
thermistor 2 100 millivolts/millimeter 
sound velocity 50 millivolts/millimeter 
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Run 5 



salinity 50 
thermistor 1 100 
thermistor 2 50 
sound velocity 20 

Rim 6 

salinity 50 
thermistor 1 100 
sound velocity 50 



Run 8 

salinity 100 

thermistor 1 100 

sound velocity 100 

Run 9 



salinity 100 

thermistor 1 100 

Run 10 

salinity 50 

thermistor 1 100 



millivolts/millimeter 
millivolts /mil lime ter 
millivolts/millimeter 
millivolts/millimeter 



millivolts /mil lime ter 
millivolts/millimeter 
millivolts/millimeter 



millivolts/millimeter 

millivolts/millimeter 

millivolts/millimeter 



millivolts/millimeter 
mi Hi volts /millimeter 



millivolts/millimeter 
mi lli volts/millimeter 



In addition to Bissett-Berman salinity measurements, 
Nansen casts were taken during each run in order to monitor 
the calibration correction which had to be applied to the 
dial readout of the base value of salinity (see Table III) . 
The correction factor was (-) 5.1 ppt. While this was a 
high correction factor, pre- experiment checks on the STD 
led to the assumption that the system was operating correctly 
in the measurement of fluctuations and was consequently used 
for that purpose. 
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Run 


Depth 

(Meters) 


Time 


Date 


Salinity 

(%) 


2 


' 4.2 


1530 


21 Oct 


33.520 


3 


9.0 


1636 


21 Oct 


33.512 


6 


4.3 


0357 


22 Oct 


33.506 


6 


4.3 


0410 


22 Oct 


33.478 


6 


4.3 


0425 


22 Oct 


33.489 


6 


4.3 


0442 


22 Oct 


33.484 . 


7 


4.3 


0552 


22 Oct 


33 . 506 


7 


4.3 


0559 


22 Oct 


33 . 480 


8 


9.5 


0654 


22 Oct 


33.459 


8 


9.5 


0659 


22 Oct 


33.438 


8 


9.5 


0715 


22 Oct 


33.471 


9 


14.6 


0725 


22 Oct 


33.469 


10 


7.8 


0810 


22 Oct 


33.455 


10 


7.8 


0817 


22 Oct 


33.463 


11 


8.2 


0839 


22 Oct 


33.453 


11 


8.2 


0852 


22 Oct 


33.455 


12 


5.7 


0934 


22 Oct 


33.440 



TABLE III 
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FIGURE 38. 



Power Spectral Density 
Thermistor 1 Run 3 
9.0m 1616 21 Oct 
Spectrum ((°C) 2 /HZ) 

Versus Frequency (0-2.5 HZ). 
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' Power Spectral Density 
Thermistor 2 Run 3 
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FIGURE 41. 

Power Spectral Density 
Thermistor 2 Run 4 
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Power Spectral Density 
Salinity Run 4 
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FIGURE 45. 

Power Spectral Density 
Thermistor 2 Run 5 
6.9m 1728 21 Oct 
Spectrum ( (°C) 2 /HZ) 

Versus Frequency (0-2.5 HZ). 
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Power Spectral Density 
Salinity Run 5 
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